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Abstract 



Data collected with Polrad swept frequency analyzer during almost 2.5 years 

of Intcrball-2 mission arc nscd to investigate the directivity properties of AKR 
observed in the Northern Hemisphere. In contrast with the previous, mostly 
statistical studies, single dynamic spectra have been analyzed. Simple modeling 
of AKR source visibiUty from the location of the observer makes possible to find 
approximate locations of remotely observed radiation sources and to model the 
overall form of the dynamic spectrum with a single set of geometric parameters 
valid during up to tens of minutes. Cases of filled and partially filled emission 
cones are found as well as smaller emission cones embedded in the main cone. 
Anisotropic AKR beams aligned with subcones as thin as 10° and sometimes 
narrower have been detected. 



1 Introduction 



Auroral Kilometric radiation has been observed since more than 30 years [Benedik- 
tov et al, 1965, Dunckel et al.„ 1970, Gurnett, 1974]. Although the body of 
available information concerning AKR is impressive there are still some very 
simple questions without a satisfactory answer. A set of basic facts, some es- 
tablished quite a long time ago, will serve as a basis for the present paper. 

1. AKR sources are located on magnetic field lines having their footprints 

within the auroral oval [Gallagher and Gurnett, 1979, Benson and Calvert, 
1979, Huff et al., 1988 and recently Gurnett et al, 2001, Schreiber et al, 
2002, Mutel et al, 2003, Mutel et al, 2004] 

2. Radiation is observed within delimited region [Green Gurnett and Shawhan, 
1977, Gallagher and Gurnett, 1979, Hashimoto, 1984]. R-X mode as the 
main product of electron maser cyclotron mechanism commonly accepted 
for AKR undergoes strong upward refraction close to the source [Calvert, 
1981, de Feraudy and Schreiber, 1995, Louarn et al., 1996a]. 

3. AKR is confined to the seemingly filled [Green et Gallagher., 1985, Green 
et al, 2004] emission cone, on the other hand cyclotron maser theory 
predicts strongest amplification of the radiation perpendicularly to the 
local magnetic field in the source [Wu and Lee, 1979, Pritchett et al., 
2002] what can suggest a hollow emission cone. 

4. AKR dynamic spectra contain a lot of features, some arc rather long-lived 
(tens of minutes) and aligned with the lower boundary of the dynamic 
spectrum 

1.1 Problems to be addressed 

Present paper focuses on four primary questions: 

1. What part of AKR radiating auroral oval can be seen from the single 

spacecraft location? 

2. How to determine approximate AKR source region location using the 
shape of dynamic spectrum? 

3. Is the AKR omission cone filled, hollow (or both)? 

4. What can be inferred on the AKR beaming from the dynamic spectrum? 

Answers to all above questions rely on simple geometrical reasoning and as it 
will be shown later, the lack of detailed information about the remote AKR 
sources and their exact positions will not affect the general conclusions of the 
paper. 

2 Instrumentation and data 

AKR dynamic spectra analyzed in the paper belong to the data set collected dur- 
ing almost 2.5 years of Interball-2 mission effective life. Satellite was launched 
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on August 29, 1996 into an elliptical orbit with an inclination of 62.8°, an apogee 

altitude of 19140 km, and a perigee altitude of 772 km. Dynamic spectra were 
recorded by step frequency analyzer Polrad [Hanasz et al, 1998] swept over one 
of two frequency ranges: 4 kHz - 0.5 MHz, or 4 kHz - 1 MHz, with repetition 
periods of 6 or 12 s and a frequency resohition of 4 kHz. Three independent 
channels allowed reception of signals from one or three perpendicular electric 
antennas {X - monopole 11 m long, Y and Z - dipoles 22 m long). Antenna 
system was spinning around the axis parallel to the X antenna with a period 
of 120 s. Polrad was working in one of its operational modes as a polarimeter. 
One frequency step corresponds to 25 or 50 ms collecting time with 6 ms output 
filter time constant. The orbit was chosen in such a way that the satellite path 
in its apogee part was roughly tangent to the auroral oval. Data were recorded 
during northern passes when spacecraft was above Earth's radiation belts. 

3 Choice of events 

Dynamic spectra with contiguous emissions lasting for tens of minutes were se- 
lected. lnterball-2 was not crossing the AKR source regions but spectra recorded 
during source approach and subsequent receding are readily available. Most in- 
teresting are spectra with aligned bottom and upper frequency edges. It will 
be shown that they contain information about filling of the emission cone. On 
the other hand, bright features lasting for some time (at least a few minutes) 
can carry information about AKR beaming. Necessary care was taken to avoid 
proximity of plasmasphere. Refraction and/or reflection caused by plasmasphere 
obscuring AKR view has not been taken into account in proposed model. 

4 Model 

Two scenarios of AKR source observations are employed (in both cases R-X 
mode has been assumed as a basic carrier of AKR energy). 

1. For calculating AKR visibility from a given point of the Interball-2 orbit 
a set of radiating magnetic field lines anchored to the auroral oval defined 
by its invariant latitude is used. Field lines fill the whole 24 hours MLT 
range. A rectilinear AKR propagation between source and observer is 
assumed. 

2. For finding AKR source position on the auroral oval and for emission cone 
investigation the source is located on a single magnetic field line. Both 
the rectilinear propagation to the observer and a simple raytracing (based 
on Calvert's, [1981]) formulas have been taken into account. 

Earth's magnetic field is approximated by a centered dipole. Positions of the 
AKR sources on individual magnetic field lines are defined assuming escape of 
EM radiation slightly above the local R-X mode cutoff frequency fx [Wu and 
Lee, 1979, de Feraudy and Schreiber, 1995, Louarn et Le Quean, 1996a, Ergun 
et al., 2000]. The ratio of AKR escape frequency / to fx is a free parameter 
of the model. Second parameter relates to the constant A^o in Nsumei et al., 
[2003] polar electron density model adopted for this paper. Maximum AKR 
emission frequency range 0.025 - 0.800 MHz has been adopted, practically it 
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was limited to the observed AKR bandwidth. Frequency range was divided into 

250 frequency steps, corresponding points on magnetic field lines were treated 
as nodes with attached AKR emission cones. Cone axis was tangent to the 
magnetic field line. 

4.1 AKR source visibility from a given satellite position 

Simple "visibility" diagram for the AKR is proposed. Assuming that the AKR 
can be produced within the whole MLT range of the auroral oval (in this sim- 
plified model labeled by a single invariant latitude value) one can construct 
for a given position of the observer a contour plot of the propagation angle a 
(between direction to the observer and the direction antiparallel to the local 
magnetic field in the source B) as a function of emission frequency (or vertical 
position of the source on the field line) and MLT of that line (Fig.l). Having 
in mind that R-X mode in the source region is generated most efficiently at 
angles « 90° to B [Hilgers, de Feraudy and Le Queau, 1992, Ergun et al, 1998, 
Ergun et al., 2000], and refracted upward in the vicinity of the source, radiation 
cannot move out of a half-sphere. In the real situation a < 90° [Bahnsen et al., 
1987, Roux et al., 1993, Louarn et Le Queau, 1996a]. R-X mode must climb up 
the generation cavity before escaping to the free space. That leads to upward 
rotation of k vector. Similar additional effect can be expected as a result of 
refraction at the cavity edges. Therefore line a = 90° delimits on / — MLT 
plane region of AKR visibility. 

The trivial property of the diagram is its symmetry with respect to the 
observers position in MLT. It results from assumption of constant invariant 
latitude for all parts of the oval. Another property: the lowest observable 
frequency cannot be lower than the R-X mode cutoff at the observers location. 
For remote observations it is always higher than that. Partially filled emission 
cone, when AKR is radiated for a range of emission angles, will map for fixed 
source position to the limited frequency range on the diagram. For a given 
frequency a contiguous part of the oval (or two segments in case of partially 
filled emission cone) can eventually host the AKR source. 

"Visibility" diagram can be used: 

1. for single satellite observations to decide what parts of the oval eventually 
radiating AKR can (or cannot) be seen from the observer's location. Such 
diagrams can be quite distinct for different satellites: low orbiting FAST 
will "see" at a given position only a few hours segment while POLAR could 
gather AKR radiation generated probably in many sources over most of 
the oval 

2. for simultaneous observations from several spacecraft to define intersec- 
tion of visibility diagrams and locate regions on the oval producing AKR 
eventually accessible for all spacecraft. Different satellites are usually hit 
by different beams of AKR and the intersection of their visibility diagrams 
may correspond to different propagation angles and azimuths 

3. for approximate determination of emission cone opening angles (taking 
into account the instantaneous AKR frequency range). 

Keeping in mind simplicity of the proposed approach such results should be 
treated as informative but sometimes they can be quite useful, for example for 
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restricting the positions of possible AKR source regions or for taking decision 
whether remotely observed AKR is produced on dayside or nightside part of the 
oval. 

4.2 Modeling of the shape of AKR dynamic spectrum and 
finding position of the source on auroral oval 

AKR "visibility" diagram offers a simple geometrical information concerning 
the radiation sources when there is no good data about the source position on 
the oval. Having (or assuming) that position it is possible to reproduce overall 
shape of AKR dynamic spectrum. Every modeled emission cone consists of a set 
of embedded subcones parametrized by their propagation angle a. For upgoing 
rays a is smaller than 90°. The basic geometric outline is demonstrated in 
Fig. 2. For the raytracing purpose angle d (between k vector and the direction 
antiparallel to the local magnetic field in the source B) analogical to a will be 
introduced. 

For equally spaced points on the orbit (in this paper one minute apart) 
and frequencies (250 steps for dynamic spectrum frequency range) one has to 
calculate for the rays reaching the satellite the values of a and Q angles at 
the exit of the radiation from AKR source. Angle a is calculated from the 
dot product of the B and direction to the observer. More complicated task 
of finding d angles is based on the interpolation between many rays launched 
from the source. Curves of constant a (and &) angle will be overlaid on the 
dynamic spectrum. Fig. 3 shows outline of projection of two a — const curves 
onto the dynamical spectrum. In the present paper analytic raytracing has been 
based on formulas published by Calvert [1981]. His dispersion relation for the 
extraordinary waves in a low-density plasma has the form: 

nl + {R+l)-nl/2 = R (1) 

where 

R=l-X/{l-Y) 

Y = fce/f 

and fpe is the plasma frequency, fee equals electron gyrofrequency and / is 
the AKR frequency. 

Refraction index components and are correspondingly perpendicular 
and parallel to the magnetic field. For rays launched perpendicularly to the 
magnetic field Uz — at the beginning of ray path. Equation(^) defines Rq 
value at that level: 

Ro^iRo + 1)-^ (2) 

More general form of equation should explicitly take into account oblique 
propagation with nonzero ti^; and starting components. 
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Solving ||2J) with respect to and substituting the result to ^ gives: 




i?0 



(4) 



As in Calvert's paper raytracing is based on the Snell's law and the differen- 
tial equation for ray paths will result from the equality dz/dx = Uzjux- There 
are two consequences of this approach: 

1. no lateral deviations of the rays are taken into account 

2. wave normal and group velocity vectors are supposed to be parallel, this is 
certainly true close to the directions parallel or perpendicular to the mag- 
netic field, for the intermediate angles and the set of plasma parameters 
adopted in the paper the angles difference can attain up to a few degrees. 
As in Calvert's paper comparisons were made with results of 3-D raytrac- 
ing code based on Haselgrove set of differential equations \de Feraudy and 
Schreiber, 1995], and no significant discrepancies were found. 

Using Calvert's substitutions for X and Y in R definition {X — const, 
Y ^ 1 — A - d and z — A ■ d/X where A is magnetic field gradient and d distance 
along field line) equation (0J can be reformulated to its final form: 



where zq = — i?o) 

This differential equation has a rather cumbersome closed form integral sub- 
sequently used for ray paths calculation. 

Rays span the full 360° azimuthal angular range around radiating magnetic 
field line. For every ray arriving at the satellite its azimuth with respect to the 
magnetic shell surface element has been determined (starting from West through 
South to East). A set of constant azimuth values curves has been overlaid 
on the dynamic spectrum, results are shown in Fig. 4 and are discussed later. 
There can be problems with determination of absolute azimuth in case of folded 
source region (see scenario proposed by Louarn and Le Queau, 1996a, 1996b or 
Pritchett et al., 2002) but the azimuth values together with the range of a give an 
idea about extension of the AKR beam producing observed dynamic spectrum. 
Refraction plays an important role in rays deflection close to the AKR source [de 
Feraudy and Schreiber, 1995] and in general case of the spacecraft not crossing 
the source region the 9 curves will be placed on the dynamic spectrum above 
the a set. The a curves can be subsequently used for the estimate of emission 
cone opening angle. It is possible to look for eventual dependency of emission 
cone opening angle on the frequency (within the frequency band of the dynamic 
spectrum). Fits of the curves to the shape of the dynamic spectrum, especially 
to its lower frequency boundary, depend in this model on four parameters: Nq 
and f/fx, the invariant latitude and MLT of the source region. Initial guess 
of source MLT coordinate relates to satellite MLT for the lowest point of the 
spectrum (if bottom part of the spectrum has a typical saucer-like form). A^o 
controls mainly vertical position of the curves and to some extent curvature. 



dx 
dz 



2 • z • (zo - 1) 



(5) 



z + 2- Zq - {ta.n9) ^ • (z - 1) - zq 
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In the present model A^o varies within the Umits 500 - 2000 cm~^ and f/fx 
between 1.01 and 1.10 (based on estimates by de Feraudy and, Schreiber, 1995 
and Louarn and Le Queau, 1996b). Changes of invariant latitude act in more 
complicated way: they can alter curvature but also rotate curves and change 
vertical position of the whole set. Similar transformations can be observed 
while changing MLT value. For cases presented in this paper good fits can be 
produced within a box ±0.5-j- 1.0 hour in MLT (depending on the case) and ±2° 
in invariant latitude. For a fixed source position on the oval variations of iVo 
and f/fx change only slightly positions of a curves while the 9 curves can move 
appreciably within the frame of the dynamic spectrum. It can be expected for 
the source region height determined mainly by the approximate equality of the 
AKR frequency and the local electron gyrofrequency fee in the source region. 
That means relatively good determination of emission cone opening angle even 
without good knowledge of the f / fx, and iVo values and of initial k vector 
direction. On the other hand, the fit of 9 curves to the dynamic spectrum takes 
into account R-X mode refraction close to the AKR source and gives additional 
refinement of the deduced source position. 

4.3 AKR emission cone: is it fully or partially filled? 

Interball-2 spacecraft trajectory is often tangential to the auroral oval. Two 
basic geometrical situations can be expected: 

1. tangent crossing of the filled emission cone 

2. tangent crossing of the empty (or partially filled) emission cone. 

Both situations are demonstrated in Fig. 3. For filled emission cone the upper 
frequency (dashed line) does not change much with time (assuming approxi- 
mately stable position of the AKR source bottom). But if the empty part of 
emission cone will be traversed, it can be noticed in the overall form of the dy- 
namic spectrum (continuous line representing upper spectrum boundary). The 
bottom and upper part of the spectrum will be aligned, there will be no constant 
upper frequency bound. Lines a = const overlaid on the dynamic spectrum can 
be used to measure actual size of the emission cone. Their intersections with 
lines / = const gives information about the size of the emission cone produced 
at a given frequency and traversed by the observer. In most cases satellite 
will cross only part of the emission cone, if this part is filled, no conclusion 
concerning the remaining part of the cone can be reached. 

4.4 Structures in the spectrum and the AKR beaming 

There are still two other interesting possibilities: 

1. narrowband structures inside the dynamic spectrum aligned with a = 
const OTc 6 = const lines 

2. narrowband structures at the bottom of the dynamic spectrum 

First case corresponds to the AKR emissions beamed into cone with narrow 
edges, extension in time can correspond to the natural duration of the process 
generating the structure or to the azimuthal extension of the beam. Structures 
at the bottom of the dynamic spectrum can represent caustics produced for the 
rays emitted close to ^ = 90° [de Feraudy and Schreiber, 1995]. 
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5 Results 



Some representative cases illustrate in Fig. 4 the classification based on the 
model. 

(a) 24/07/98 0600-0700 UT. This is a "typical" Polrad's spectrum with 
rather wide emission cone, probably partially empty and attaining ~ 35° wide 
filling. Bright (more than 20-30 dB above surrounding background) wideband 
emissions lasting for about 20 minutes and corresponding to a beam about 
10° — 20° wide in opening angle are also visible. 

(b) 15/08/97 0600-0700 UT - about 40° - 50° wide very strong structure 
at the beginning of dynamic spectrum. Azimuthal range is not well defined and 
slight changes ±0.5 hour in MLT can move the azimuthal lines pattern up to few 
tens of degrees since satellite is moving close to the AKR source region. The low 
frequency part resembles broadband electrostatic emissions, but most probable 
explanation is the nonlinear response of the receiver to the very strong signal. 
This interpretation is supported by weak harmonic structures visible above the 
most intense part of the dynamic spectrum. 

(c) 29/03/98 1700-1800 UT - quite extended emission angles range (about 
50° with empty cone between 0° and « 20°), but with a 20° azimuthal gap 
around 70°. It is worthwhile to note nice alignment of constant azimuth lines 
with vertical borders of the dynamic spectrum. The stepwise rise of the dy- 
namic spectrum lowest frequency (~ 43min) illustrates the sudden change of 
the R-X mode emission cone opening angle, the much weaker background radi- 
ation below propagates in L-0 mode (unpublished results of Polrad polarimeter 
provided by M. Panchenko) 

(d) 20/07/97 1300-1400 UT - narrow empty emission cone about 15° 
thick, bright structure inside less than 10° wide, up to four fine structures as 
thin as 1° — 2°, quite extended (about 60°) azimuth range. Dayside source. 

(e) 30/10/97 0000-0100 UT - empty emission cone about 10° thick em- 
bedded inside diffuse AKR radiation. 

(f) 27/05/97 1600-1700 UT - two separate AKR emissions containing 
structures aligned with 6 = const lines about 2° and 10° wide. 

Visibility diagrams for all above mentioned cases at the satellite closest ap- 
proach to the source are shown in Fig. 5. The shaded area corresponds to the 
frequency ranges for which AKR is not seen in the spectrum. Table 1 presents 
best fit parameters. 

6 Discussion 

Averaged intensity variations across AKR emission cone were extensively stud- 
ied by Green et Gallagher, 1985. Conclusions based on such statistical approach 
generally favor filled emission cone although with different opening angles vary- 
ing with frequency. Latest findings by Green et al., 2004 show seasonal and 
solar cycle variations of the AKR dynamic spectra additionally complicating in- 
terpretation of emission cone geometry. In such a situation statistical approach 
wipes out eventual individual properties of the AKR emissions. There were 
attempts to interpret some peculiarities of Green and Gallagher, 1985 results, 
especially for lower frequencies (56.2 kHz and 100 kHz), as pointing toward 
empty emission cone [Calvert, 1987 or two-beam emission pattern Zarka, 1988]. 
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The overall form of all dynamic spectra presented in the paper can be delin- 
eated by a single set of a. lines. Intcrball-2 missed AKR source regions and cal- 
culated positions of the sources are approximate. Positions confined to an error 
box ±1.0 hour in MLT and ±2° in invariant latitude implicate relatively small 
source regions stable during at most tens of minutes although sometimes differ- 
ent frequency bandwidths can be excited. They agree with the present estimates 
made by Mutel et al., 2004 for the short AKR bursts. They report position of 
AKR source ccntroid wandering during 15 minutes within a region constrained 
by 2100 MLT ±1 h and 70° ±3° invariant latitude. The curves a = const fitting 
the bottom limits of the dynamic spectra give cone half-opening angles for the 
dynamic spectrum frequency band. Highest a value reported by Louarn and Le 
Queau, 1996a attained about 70°. For cases presented in this paper they vary 
between almost 90° (when the satellite is located close to the source region) and 
60°. The values close to 90° may be misleading. The sets of curves overlaid on 
the dynamic spectra are calculated for the best fit of filamentary source aligned 
with the specified magnetic field line. If the sources are scattered within the 
error box then real propagation angles can be slightly different because close to 
the source region refraction may still bend rays [de Feraudy and Schreiber, 1995] 
and the estimated propagation angles will differ from those measured at greater 
distances. For remote sources the lines corresponding to 90° and 80° are always 
much closer to each other, than the lines representing smaller angles. For the 
analysis of structures within the dynamic spectrum an error in determination of 
the cone opening angle does not matter much, structure 10° wide can shrink to 
say 7.5° or expand to 15° depending on their vertical positions on the dynamic 
spectrum but this is still a relatively narrow beam. For AKR sources close to 
the satellite as in Fig. 4b and Fig. 4c both sets of lines do not differ much. 

Good alignment of a = const lines with bright structures in dynamic spec- 
trum implies constant opening angle of the emission cone (or embedded subcone) 
within the frequency band occupied by the structure. 

Azimuthal extension of the radiation beam corresponding to the bright struc- 
tures visible in analysed AKR spectra vary from a few degrees to more than 80°. 
At the present stage it is not possible to separate space and time effects and in- 
dependent observations from other spacecraft located in that same beam would 
be helpful. But at least the minimum azimuthal extension of the beam can 
be estimated. It is difficult to compare our beam sizes with AKR short bursts 
beams discussed recently by Mutel et al, 2004 because our findings suggest 
rather small ~ 10° latitudinal thickness of the bright parts of emission cone 
what suggests anisotropic beams. Mutel et al., 2004 comments on AKR bmst 
beams sizes in the frame of Cluster spacecraft separation as rarely (10 ^ 30%) 
reaching 20° but apparently does not take into account possible selection effects 
introduced by orientation of the baseline projection on the emission cone edge. 

Some fine structures {'^ 1° — 2° wide) observed by Polrad can be interpreted 
as produced by thin conical beams located on the wall of the hollow emission 
cone. They resemble thin Jovian DAM arcs recently discussed by Kaiser et al, 
2000. 
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7 Conclusions 



Dynamic spectra investigated on case by case basis contain a lot of information 
concerning radiation diagram of AKR complementary to information gathered 
in previous statistical studies. Simple model presented above is sufficiently good 
to catch and reproduce important characteristics concerning AKR directivity. 

1. General shape of the AKR dynamic spectrum can be modeled during tens 
of minutes with one set of model parameters what implies stability of the 
overall source size and position within the ±1 h and ±3° error box during 
that period of time. It is possible to determine approximate position of 
the AKR source region even for remote observations. 

2. In many situations AKR emission cone is not filled, judging from the cases 
discussed in the paper, "walls" can be 10 -j- 60° thick. 

3. Some strong (20 db and more above the AKR background level) and 
relatively narrowband (a few tens of kHz wide) structures lasting for about 
ten minutes correspond to thin parts of emission cone (» 10° wide). Other 
structures a few kHz wide can be produced by very narrow beams attaining 
1° — 2° like those reported for Jupiter DAM arcs. 

4. Structures in the dynamic spectra aligned with constant emission cone 
opening angles are consistent with the presence of anisotropic AKR beams 
(subcones) restricted both in elevation and in azimuth. 

5. Alignment of the curves of constant propagation (a) and emission {9) an- 
gle with low frequency borders of the dynamics spectrum or with bright 
structures inside the spectrum implies constant opening angle of the emis- 
sion cone for the range of frequencies covered or corresponding altitude 

range of the source regions. 

6. There are dynamic spectra showing gaps in time what can eventually be 
interpreted as a result of emission cone "bitten out" in azimuth. Confir- 
mation of this conjecture needs however multispacecraft observations. 
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Figure 1: Example of AKR visibility diagram for Polrad on 24/07/98 0630 UT. 
Hatched regions correspond to the parts of the auroral oval not visible from the 
momentary observer location 



Figure 2: AKR emission cone/embedded subcone concept 
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27/05/97 


70° 
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1.10 


500 CTO"^ 



Table 1: Source peirameters 



Figure 3: AKR emission conc/subcones and beams aligned with cone edges 
projected on the dynamic spectrum. 
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Figure 4: AKR dynamic spectra with superposed = const - solid line, a = 
const - dashed line and azimuth angle - dotted line 



Figure 5: Visibility diagrams for the spectra shown in Fig. 4 at the satellite 
closest approach to the source. The shaded areas correspond to the frequencies 
not seen in the spectrum, hatched regions cannot be seen at the given observer 
location 
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